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Intricate dynamics of a deterministic walk confined in a strip 
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Abstract. - We study the dynamics of a deterministic walk confined in a narrow two-dimensional 
space randomly filled with point-like targets. At each step, the walker visits the nearest target 
not previously visited. Complex dynamics is observed at some intermediate values of the domain 
width, when, while drifting, the walk performs long intermittent backward excursions. As the 
width is increased, evidence of a transition from ballistic motion to a weakly non-ergodic regime is 
shown, characterized by sudden inversions of the drift velocity with a probability slowly decaying 
with time, as 1/t at leading order. Excursion durations, first-passage times and the dynamics of 
unvisited targets follow power-law distributions. For parameter values below this scaling regime, 
precursory patterns in the form of "wild" outliers are observed, in close relation with the presence 
of log-oscillations in the probability distributions. We discuss the connections between this model 
and several evolving biological systems. 
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Introduction. — Deterministic walks in disordered 
environments have received an increasing attention over 
the past years. They describe diffusion processes follow- 
ing non-random rules and have applications, among oth- 
ers, to the study of the displacements of individuals in 
complex landscapes. Examples are human travels [1,2], 
human displacements in a city [3], movement patterns of 
hunter-gatherer [4] or foraging animals [5] . 

From a given position, the next site visited by a purely 
deterministic walker is assigned from a given set of rules 
and not stochastically. These walks still have probabilistic 
and fluctuating features if the environment is random or 
heterogeneous. Interesting dynamics have been observed, 
such as normal [6] or anomalous diffusion [7], behaviors 
analogous to that of the Lorentz gas [8] , cycles with power- 
law distributed periods [1,9] or Levy-like step length dis- 
tributions [10,11]. Complex behavior can emerge from 
very simple rules, e.g, when each individual step optimizes 
a given cost function. Some properties of deterministic 
walks have also been used as tools to process large data 
sets in galaxy surveys [12], thesaurus graphs [13] or for 
pattern recognition [14]. 

In many situations, in particular biological, the deter- 
ministic walker itself changes the medium, which intro- 
duces memory [8]. An important case is the self-avoiding 
walk (SAW), which can be implemented to model biolog- 
ical systems with negative feedbacks that tend to avoid 



past behaviors. A simple example is that of a foraging an- 
imal relying on mental maps to navigate an environment 
composed of food patches that are not revisited after they 
have been depleted [5]. In a different context, the brain 
activity has been modeled by random walks keeping mem- 
ory of their complete history in order to avoid persistent 
patterns; recent memory loss producing pathological rep- 
etitions, like in the Alzheimer's disease [15]. In evolution- 
ary ecology, the well-known Red-Queen principle assumes 
that any organism must constantly evolve in order to pre- 
vent its predators or preys to adapt to an otherwise pre- 
dictable behavior. Similar considerations can apply to the 
dynamics of technological innovations [16]. Freund and 
Grassberger introduced some time ago a self-avoiding de- 
terministic walk model in disordered two-dimensional do- 
mains, mimicking evolving organisms in phenotype land- 
scapes [7] . These kinds of models are very difficult to han- 
dle analytically; they are firstly dynamical and usually not 
equivalent to canonical SAWs [17]. 

Here, we study a minimal model of a deterministic walk 
with a SAW constraint (in the infinite memory limit) and 
confined in a nearly one dimensional random medium. In 
an evolutionary context, whereas SAWs can be justified by 
natural selection, organisms also have developmental con- 
straints due to limited phenotypic variability [18]. Simi- 
larly, the development of human artifacts is restricted by 
design limits. We model this important constraint by a 



D. Boyer 



narrow random medium where the walker can evolve with- 
out bounds only in one direction. In narrow landscapes, 
the model exhibits very rich dynamical features not ob- 
served in unbounded ones, such as intermittent behavior, 
scaling laws, discrete scale invariance and very large events 
(outliers). A discussion of these results is then presented. 

Model description. — Consider a two-dimensional 
strip of width I and infinite length along the horizontal 
direction. The strip is randomly filled with fixed point- 
like targets with uniform number density po, representing, 
say, food patches for a foraging animal or phenotypes for 
an evolving species. The only control parameter is the re- 
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duced domain width, defined as 8 = l/lo, with Iq — p Q 
the characteristic distance between neighboring targets. 
At time t = 0, a walker is located at some target with 
coordinates (xq, yo), taken as the origin. Two rules of mo- 
tion are then recursively applied: the walker (i) moves to 
the nearest available target, (ii) does not visit a previ- 
ously visited target. When the new target is reached, t is 
updated to t + 1. 

The medium can be made one(two)-dimensional in the 
limit i « 1 (J » 1), respectively. We will focus here 
on values of 8 of 0(1), typically in the range (2, 5), such 
that the walker has a some vertical degree of freedom but 
a practically one-dimensional motion on large scales, de- 
scribed by its horizontal coordinate x(t). In the simula- 
tions, the medium is a rectangle of area unity containing 
N targets and of width I = S/^N. Each run start near 
the middle of the domain and is stopped before the walker 
reaches the lateral vertical walls. 

Trajectories. — In the one dimensional case (8 <C 1), 
the targets are randomly distributed on a line and the 
motion is simply ballistic. After a possible short tran- 
sient, the walker breaks the right-left symmetry and al- 
ways moves to the nearest target to its right (or left) so 
that x(t) is a sum of same-sign independent random vari- 
ables with Poisson distribution. The 2c? case (8 3> 1) is 
sometime called the "tourist walk" [1,2]: the trajectories 
are not very different from 2c? random walks, although 
slightly superdiffusive [19]. 

For the cases 8 = 0(1) of interest here, the situation 
is quite different and trajectories exhibit a rich structure. 
As shown in fig. [I] at 8 = 4.1, the motion is on aver- 
age ballistic due to the confining effect of the horizontal 
walls. The numerically calculated root-mean-square dis- 
placement [xit) 2 ) 1 / 2 follows a linear behavior with time 
(not shown here). Note that the walker horizontal veloc- 
ity x(t) — x(t — 1) often changes sign: the walker per- 
forms many "backward excursions" while drifting along 
the strip. These excursions, that were observed in a pre- 
liminary study of the model [10], can be explained qual- 
itatively. A walker drifting, say, toward the left does not 
necessarily visit all the targets of a given neighborhood 
on its way and may ignore some targets. From time to 
time, rules (i) and (ii) make the walker turn back and 
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Fig. 1: Space-time diagram of a trajectory with 8 — 4.1 (x is 
in unit of Iq). 



visit these unvisited targets toward the right, until it ends 
up in a region depleted of available targets. In that case, 
a single step can brings the walker back to the unexplored 
region located to the left. 

Unexpectedly, backward excursions of all sizes can be 
observed in fig. QJi. Whereas most excursions are short, 
some can be of order 10 2 — 10 3 x ^o, e.g. near t = 9300. A 
close up of fig. (inset b) reveals further details and sug- 
gests that the trajectory is fractal. In [10,11], it was found 
that these intermittent backward excursions can lead to 
"Levy-like" distributions for the distance separating suc- 
cessively visited targets, of the form with fi ~ 1 
at 5 = 4 [11]. 

Additionally, the sign of the drift velocity can change 
suddenly at large t (e.g. at t ss 600 and t w 5000 in 
n g- Ht 1 )- Such inversions happen during a backward ex- 
cursion, at some point when the closest unvisited target 
is located, say, to the right of xq for a trajectory that was 
previously drifting toward the left. Obviously, inversions 
can not occur in the Id ballistic limit of the model. We 
investigate below the possible existence of a transition be- 
tween different dynamical regimes as 8 is varied. 

Inversion probability and first-passage times. 

— We define the explored interval at time t as 
[x m i n (t),x max (t)], where x mm (t) (x max (t)) is the coordi- 
nate of the leftmost (rightmost) visited target after t steps, 
respectively. A inversion (say, from right to left) occurs 
during the t th step if x m i n (t) — x m i n (t — 1) < and if there 
exists a time t' < t such that x max (t') — x 
and such that i mjn (f ) - x min (t n - 1) 
x max (t" - 1) = for t'<t" <t. 

We then define Pi nv (t) as the probability that an inver- 
sion (to the left or right) occurs during the t th step, and 
Po(t) as the probability that the walker crosses xq during 
the t th step. In the random walk language, Pq is analo- 
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Fig. 2: Probability that the drift velocity changes its sign at 
time t, as a function of t and for various strip widths. The 
probabilities are calculated from 8 fO 4 independent random 
media. 

gous to the probability of presence at the origin. If motion 
is essentially ballistic between two inversions, a trajectory 
crosses its origin at large times only during an inversion: 
Pinv and Po have the same asymptotic behavior. 

In Figure [H the numerically computed Pi nv (t) decays 
very slowly with time for several 0(1 ) values of the strip 
width 5. In domains as narrow as 5 = 1.5, after an initial 
steep decay, Pmv{t) exhibits a surprising fat tail. At the 
larger value S — 4.1, Pi nv (t) can be well fitted by the 
simple inverse power-law c/t, with c a constant. The same 
curves are reploted in Figure [3] as tPi nv (t) versus t: for 
5 = 4.1, the curve remains remarkably constant during 
almost 6 decades, while strong corrections to scaling are 
present below and above that parameter value (5 — 3 and 
5). Very similar results are obtained for Po(t). 

Despite that the walker crosses less frequently the origin 
than a Id random walker (where Po(t) ~ t^ 1 / 2 ), the return 
probability at large times in narrow strips (8 <C t) remains 
very high instead of being exponentially small as for usual 
ballistic motion (e.g., a Id random walker with a bias). 

Contrary to random walks, sign changes in x(t) are 
abrupt and not strongly correlated to the evolution of x(t) 
during the preceding steps (see Fig. Q}. It is therefore use- 
ful to make a connection between this result and a simpler 
two-state stochastic problem consisting of a walker moving 
ballistically on a line with two possible velocities, v and 
—v. Starting in one state, the walker change its velocity 
in the time interval [t, t + 1] with probability Pi nv (t), that 
is given. It is well known that if pi nv (t) decays faster than 
1/t, there is a finite probability that the walker remains 
indefinitely in a same state (v or —v) after reaching this 
state. If Pi nv (t) decays as 1/t or slower, the probability 
that the walker remains in a same state forever is zero. 
The behavior of the system is non-ergodic in the former 
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Fig. 3: Same data as in fig.[2j reploted as tPi nv (t) vs. t. Inset: 
details of tP inv (t) for 5 = 4 (o), 4.f (A), 4.2 (•). 

case, as the left-right symmetry is asymptotically broken, 
while it is weakly non-ergodic or ergodic in the latter 
case. Obviously, the above two-state problem only pro- 
vides an approximate description of our model and makes 
sense only in the regime where inversions are abrupt. 

We investigate more in details the possibility of a non- 
ergodic/weakly non-ergodic transition as the strip width 5 
is increased across some critical value 5 C where Pi nv (t) — 
c/t. Making an analogy between the behavior of Pinv (or 
Po) and that of a correlation function near a critical point, 
for 8 slightly below S c one may look for a standard scaling 
form: Pi nv (t) ~ t~ 1 g(t/T(S)), with g(x) a scaling function 
rapidly decaying to zero at large x and t(S) a diverging 
timescale as 5 — > S c . The inset of Fig. [3] displays tPi nv (t) 
vs. t for different values of 5 near 4.f and shows that the 
above ansatz does not hold. 

Interestingly, the probability exhibits an unusual behav- 
ior instead. First, the different curves can not be rescaled 
onto a single curve. Second, a pure power-law behavior 
was never obtained for Pi nv (t) (nor Po(t)) for the values of 
S considered in this study. Intricate corrections to scaling 
in the form of logarithmic oscillations are observed. Log- 
oscillations have been observed in a variety of systems and 
are a manifestation of the phenomenon of discrete scale in- 
variance [21]. The log-oscillations have a large period, of 
order 2 In 10, which complicates the observation of several 
periods: we can not conclude whether they converge to- 
ward a finite amplitude or are amplified. However, the 
amplitude of the oscillations is minimum at S c ~ 4.1. 

The leading 1/t decay of Pi nv at S c is probably not a 
coincidence. From a renormalization group (RG) perspec- 
tive, the model has a trivial attracting fixed point, S* = 0, 
corresponding to simple ballistic motion in very narrow 

1 In this context, weakly non-ergodic means that both states al- 
ways remain accessible to the walker, although the time interval 
between two visits diverges asymptotically [20]. 
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Fig. 4: a) Probability distribution Pi(t) of the first-passage 
time at <5 = 4.1 (A). Opposite side: probability distribution 
P(At) of the time intervals between two consecutive zeros of 
N u (t) for 5 = 4.1 (o) and S = 2.45 (•). Pi(t) and P(At) are 
obtained from 8 10 4 and 10 3 independent runs respectively, 
b) Probability distribution of the excursion duration (r > 2). 
Lines are guides to the eye. 

strips. The 1/t law indicates that S c should lie at the 
boundary of the basin of attraction of that fixed point. 
Besides, the increasing corrections from power-law behav- 
ior for S slightly above 8 C (see Fig. |3]at 5 = 5) suggest that 
RG trajectories above S c flow towards an other attracting 
fixed point (that could be S* — oo). This argument sup- 
ports the idea that the transition is not a cross-over and 
that 5 C might be a non trivial repelling fixed point. 

An other possibility is that the walk may become 
asymptotically ballistic without inversions after extremely 
large times, unreachable with standard numerical meth- 
ods. In this case, the results above would describe a very 
long transient preceding an asymptotic regime of limited 
practical relevance. 

Other insights into inversion processes can be gained 
from the distribution P\{t) of first-passage times. The 
first-passage time is defined here as the step number when 
the walker crosses for the first time xq . As shown in figure 
UK, a t" Q law with a ~ 3/2 holds remarkably well over 
nearly 6 decades in the vicinity of S c . As for Pi nv (and 
Po), log-oscillations were detected in P\ at (and near) 8 C . 

This exponent value can be qualitatively explained with 
the help of the simple two-state approximation described 
above, where the probability that the velocity changes 
its sign for the first time at time t reads: Pi(t) = 



p inv (t) exp[- JgPmv{t')dt']. If p mv ~ c/t at large times, 
then pi(t) ~ t~ a with a = 1 + c. The numerical value of 
c calculated from Pi nv at 5 = 4.1 yields a ~ 1.33. This 
value is close to, but lower than the observed 3/2. There- 
fore, inversion events are not independent but probably 
long-range correlated. It is actually surprising (and most 
likely coincidental) that the first-passage exponent is close 
to the simple value of the Id random walk [22]. 

Backward excursions and unvisited sites. We 

now come back to the description of backward excursions, 
that are much more frequent than inversion events. The 
probability distribution of excursion durations, P e (r), can 
be obtained from the sizes of the time intervals when 
(t) remains constant. As shown in figure 
in the vicinity of 5 C this distribution is also well fitted 
by a power-law behavior, P e (r) ~ t -/3 , with (3 ~ 2.7. This 
distribution has finite first moment but infinite variance. 
On average, the walker remains "trapped" in an excur- 
sion during a finite number of steps, but its progression 
is quite intermittent. For strip widths well below 5 C , ex- 
cursions are still observed and P e (r) remains fairly broad, 
although it can no longer be fitted with a power-law. Gen- 
erally speaking, backward excursions tend to restore the 
right-left symmetry of the system. For this reason they 
are reminiscent of the effect of thermal fluctuations on a 
broken symmetry phase in equilibrium. 

An other quantity of interest related to excursion dy- 
namics is the number of unvisited sites in the explored 
interval [x m i n (t), x ma x(t)}, denoted as N u (t). As shown in 
figure EH, for values of S below S c , N u (t) displays cycles of 
irregular durations analogous to oscillations in excitable 
systems. The cycles are composed of (i) a slowly increas- 
ing part on average and (ii) a fast decay down to zero. 
This behavior reflects the fact that a small fraction of sites 
are left as unvisited while the walker is drifting in the dis- 
ordered medium, leading to an increase in N u (t). These 
sites can be visited later, in a long backward excursion, 
leading to an "avalanche- type" relaxation of N u (t). (Note 
that many smaller excursions also occur during the ascend- 
ing part of N u (t).) The distribution of the time intervals 
At between two successive zeros of N u (t) is displayed in 
figure 2^. One expects P{At) and the first-passage time 
distribution P\(t) to have the same asymptotic behavior, 
as observed. At 8 — 4.1, one finds P(At) ~ (At)" 3 / 2 , 
implying that (At) = oo: N u (t) grows asymptotically un- 
bounded. In the transition region, the walker is therefore 
unable to visit all the targets of the explored interval at 
large time (it is "inefficient"). Below S c , the distribution 
P{At) decays faster (figure^) and the evolution of N u (t) 
seems to have a characteristic cut-off period. In figure 
Ek, where 5 is well below 5 C , this characteristic time is still 
very long (~ 10 ). Large avalanches, where N u drops from 
about 600 to 0, are present. 

The walker is a priori "efficient" below 5 C , since it reg- 
ularly leaves no sites as unvisited (fig. EJa,). However, 
this behavior is not persistent on very large time-scales, 



Intricate dynamics of a deterministic walk confined in a strip 




60000 



Fig. 5: Number of unvisited targets N u (t) as a function of t, 
at early (a) and late (b) times, for a same trajectory below the 
transition region (8 — 2.45). 



as shown in figure [5Jd. Surprisingly, at a given time that 
can be of order 10 6 or more, N u (t) does not come back 
to zero and starts to oscillate above a finite value. This 
happens when a large excursion fails to visit some of the 
unvisited targets left behind. After such an incomplete 
excursion, (At) obviously starts to grow with time. This 
behavior is observed in a whole range of parameter values 
below S c , down to about 5 ~ 1.30. The walk is therefore 
efficient during a finite time, until a "catastrophic" event 
with very large At occurs. 

One can draw an analogy between these "outliers" and 
very large earthquakes or financial crisis in other contexts 
[23,24]: their magnitude can not be explained from a sim- 
ple extrapolation of the distribution P(At), as displayed 
in fig. t° larger arguments. Here, outliers are unusual 
precursory patterns of the power-law distribution (in sharp 
contrast with common critical phenomena) and can be at- 
tributed to the presence of log-oscillations in P(At): As 
noticeable in figure [4^,, after an initial steeper decay, the 
distribution for S = 2.45 < 6 C very closely approaches the 
critical curve for a range of values around Ati ~ 5 10 3 , 
before decaying fast again. Although computational limi- 
tations do not allow to observe more log-periods, it is very 
likely that the next oscillation can reach (or even cross) 
the critical distribution, for some values At 2 ^> At\ (see 
figure [3] for a similar behavior at 8 = 3). Therefore, the 
emergence of very large intervals of order Ai 2 between the 
zeros of N u becomes as probable as at S c . This explains 
qualitatively the time series of figures [5^,-b, composed of 
many intervals of order 10 3 — 10 4 followed by a single one 
of much larger size, precursor of the critical regime with 
diverging (At). 

Discussion. — We have shown that a simple self- 
avoiding process taking place in a confined Poissonian ran- 
dom medium can display complex dynamics and broad 
distributions in a wide parameter range. Quenched disor- 



der introduces randomness in the model, that otherwise 
follows simple deterministic rules. Similar results as re- 
ported here should be observed in a semi-infinite strip with 
the walker initially located at one end, with the difference 
that the walker would drift without inversions. 

In ref. [7], the trajectories generated by the Red Queen 
rules can be similar to random walks after time scales that 
depend strongly on the lattice size and geometry [7]. In 
contrast with the Red Queen Walk, where sites can be 
revisited after a very long time, our model has infinite 
memory, leading to intermittent and complex behavior. 

The evolution of single species is known to be intermit- 
tent and not gradual, long period of stasis being "punc- 
tuated" by burst of rapid biological changes [25]. Such 
active periods might be driven by the internal dynamics 
of evolution. According to the fossil record, the number 
of genera with a lifetime r follow a power-law N ~ 
with (3 ~ 2 [26]. The Bak-Sneppen model [27] considers 
interacting species with high mutation barriers, leading 
to self-organized critical states with (3 — 1.1 [26]. This 
evolution is slower than observed because it occurs by col- 
lective modes, or avalanches. Changes are easier in our 
model (where explicit interactions are ignored), but still 
intermittent. The distribution of time-intervals between 
successive changes in x max , for instance, is fitted with an 
exponent j3 = 2.7. One may speculate that phenotypic 
restrictions could play a role on the punctuated dynamics 
of evolution, in addition to species interactions. 

Our system does not become critical in an ordinary 
way. At a critical width 8 C , the inversion probability 
of the drift velocity decays as a power law with small 
log-periodic corrections. In a first harmonic approxima- 
tion [21], P inv (t) ^ ct _1 [l + ai cos(27rlnt/lnA)], with 
Oi < 1 and A ~ 100. This asymptotic regime is nu- 
merically hard to reach, as observed in other problems 
with log-oscillations [15,30]. The leading term above pre- 
cisely represents the law that separates, in analogy with a 
two-state stochastic process, asymptotically ballistic (non- 
ergodic) trajectories and walks that keep changing direc- 
tion indefinitely. 

The log-oscillations present in various distribution func- 
tions indicate the presence of a hierarchy of time-scales 
related to each other by a particular scaling factor A, such 
that P(Xt) ~ X a P(t). These oscillations are often dis- 
played by cooperative phenomena taking place on hierar- 
chical structures (spin models near criticality [21], contact 
processes [28]), or by random walk models with memory 
[15,29], among other examples. In ref. [30] log-oscillations 
appear in a simple biased Id random walk model in a 
disordered medium containing a small fraction of "slow" 
sites, where the walker jumps in the direction opposite to 
the bias with a probability close to one. These slow sites 
are somehow analogous to our (dynamically generated) 
backward excursions. A crucial ingredient leading to log- 
oscillations in [30] is spatial discreteness, where clusters 
of slow sites trap the walker during a time that increases 
exponentially with the cluster size. In our model, no such 
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discreteness is apparent. Instead, inversion events are cor- 
related in a complicated way: if the strip is sufficiently 
narrow, an inversion can not occur in the interval [t±, 2t±] 
if an inversion occurred at time t\. The time intervals 
separating inversions might introduce a particular scaling 
factor, although its precise origin is unclear. 

Well below the critical region, the distributions exhibit 
log-oscillations of irregular amplitudes whose maxima can 
be identified with outliers, that are "wild" precursors of 
critical fluctuations. From the above discussion, the char- 
acteristic size of these events can be roughly extrapolated 
as being of order A™ (n = 1,2...), and as probable as at 
criticality. At <5 C , these specific scales are mixed with all 
the others (ai C 1), in a practically scale invariant distri- 
bution. A detailed study the behavior near 5 C remains to 
be done. Correlation functions (e.g. velocity) other than 
Pi nv might exhibit clearer scaling relations. 

In this scenario, it is however clear that outliers exist in 
a wide parameter range and do not even require that the 
bulk of the distribution follows a power-law. On the con- 
trary, they are off-critical events by nature. This property 
has to be contrasted with more common views in seismol- 
ogy, for instance, where outliers are either considered as 
coming from the tail of power-law distributions [31], or, in 
a more refined way, as coming from a bump at large sizes 
in a otherwise power-law distribution of bulk events [32]. 

* * * 

Fruitful discussions with H. Larralde, F. Leyvraz, 
M.G.E da Luz, O. Miramontes, A. Robledo and G.M. 
Viswanathan are gratefully acknowledged. 

REFERENCES 



[12] Elson E.C., Bassett B.A., van der Heyden K. and 
Vilakazi Z.Z., Astron. & Astrophys., 464 (2007) 1167. 

[13] Kinouchi O., Martinez A.S., Lima G.F., Lourencjo 
G.M. and Risau-Gusman S., Physica A, 315 (2002) 665. 

[14] Campiteli M.G., Batista P.D., Kinouchi O. and 
Martinez A.S., Phys. Rev. E, 74 (2006) 026703. 

[15] Cressoni J.C., da Silva M.A.A. and Viswanathan 
G.M., Phys. Rev. Lett, 98 (2007) 070603. 

[16] Kauffman S.A., At home in the universe (Oxford Uni- 
versity Press, Oxford) 1995. 

[17] Ponmurugan M., Narasimhan S.L. and Murthy 
K.P.N. , Physica A, 371 (2006) 171. 

[18] Arnols S.J., In: Behavioral Mechanisms in evolutionary 
ecology, edited by Real L.A. (University of Chicago Press, 
Chicago) 1994, p. 258. 

[19] Lopez-Corona 0. and Boyer D., unpublished. 

[20] BOUCHAUD J. P., J. Phys. I France, 2 (1992) 1705. 

[21] Sornette D., Phys. Rep., 297 (1998) 239. 

[22] Redner S., A guide to first-passage processes (Cambridge 
University Press, Cambridge) 2001. 

[23] Sornette P.. larXiv:0707.2194H . 

[24] Sornette D., In: Extreme Events in Nature and So- 
ciety, edited by Albeverio S., Jentsch V., Kantz H. 
(Springer, Berlin) 2006, p. 95. 

[25] Gould S.J. and Eldredge N., Paleobiol, 3 (1977) 115. 

[26] Sneppen K., Bak P., Flyvbjerg H. and Jensen M.H., 
Proc. Natl. Acad. Set. USA, 92 (1995) 5209. 

[27] Bak P. and Sneppen K., Phys. Rev. Lett, 71 (1993) 
4083. 

[28] Bab M.A. and Albano E.V., J. Phys. A-Math. Theor., 

41 (2008) 045001 

[29] Kenkre V.M.. IarXiv:0708.0034k 2. 

[30] Bernasconi J. and Schneider W.R., J. Phys. A: Math. 

Gen., 15 (1982) L729. 
[31] Bak P., Tang C. and Wiesenfeld K., Phys. Rev. Lett, 

59 (1987) 381. 

[32] Gil L. and Sornette D., Phys. Rev. Lett, 76 (1996) 
3991. 



[I] Lima G.F., Martinez A.S. and Kinouchi O., Phys. Rev. 
Lett, 87 (2001) 010603. 

[2] Stanley H.G. and Buldyrev S.V., Nature, 413 (2001) 
373. 

[3] Chowell G., Hyman J.M., Eubank S. and Castillo- 
Chavez C, Phys. Rev. E, 68 (2003) 066102. 

[4] Brown C.T., Liebovitch L.S. and Glendon R., Hum. 
Ecol, 35 (2007) 129. 

[5] Boyer D., Ramos-Fernandez G., Miramontes O., 
Mateos J.L., Cocho C, Larralde H., Ramos H. and 
ROJAS F., Proc. R. Soc. B, 273 (2006) 1743. 

[6] Grosfils P., Boon J. P., Cohen E.G.D. and Buni- 
movich L.A., J. Stat. Phys., 97 (1999) 575. 

[7] Freund H. and Grassberger P., Physica A, 190 (1992) 
218. 

[8] Bunimovich L.A., Physica D, 187 (2004) 20. 
[9] Derrida B. and Flyvbjerg H., J. Physique, 48 (1987) 
971. 

[10] Boyer D., Miramontes O., Ramos-Fernandez G., 
Mateos J.L. and Cocho G., Physica A, 342 (2004) 329. 

[II] Santos M.C., Boyer D., Miramontes O., 
Viswanathan G.M., Raposo E.P., Mateos J.L. 
and DA Luz M.G.E. , Phys. Rev. E, 75 (2007) 061114. 



